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1. Introduction 

In 1972 Montgomery [15] and Dyson [12] discovered that pairs of zeros of the Riemann 
zeta-function are distributed like pairs of eigenvalues of random unitary matrices. Part of 
this discovery could be proven, under the assumption of the Riemann Hypothesis, and part 

l 
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relies on a heuristic based on the Hardy-Littlewood conjectures for the distribution of prime 
pairs. 

Odlyzko [16] , in the 1980s carried out a substantial numerical test of Montgomery's con- 
jecture which provided stunning visual substantiation. 

Subsequently Rudnick and Sarnak [18] showed that the limit high on the critical line of the 
n-correlation of the zeros of the Riemann zeta-function agreed with that of unitary matrices, 
provided that the test function had a Fourier transform with limited support. 

At the same time, Bogomolny and Keating [2] showed how the Hardy-Littlewood conjec- 
tures could be used to derive the asymptotic limit of the n-correlation. 

Bogomolny and Keating [lj also investigated the difference between Montgomery's limiting 
pair-correlation conjecture and the data of Odlyzko. A close examination of Odlyzko's 
data revealed that lower order terms, likely to be of an arithmetic nature, were present. 
They derived formulae for these lower order terms, initially using the Hardy-Littlewood 
conjectures, but subsequently developing a method whose point of departure was the trace 
formula of Gutzwiller. They gave full details for the lower order terms for the 2-point 
correlation, as well as numerics showing the goodness of fit, whereas for three and higher 
correlations, they outlined several methods which lead to these lower order terms. 

In this paper, we present a different approach to obtaining these lower order terms for 
n-correlation. Our approach is based on the 'ratios conjecture' of Conrey, Farmer, and 
Zirnbauer [7J [5] (see also [H]). Assuming the ratios conjecture we prove a formula which 
explicitly gives all of the lower order terms in any order correlation. (In the final section we 
write down the first four correlations.) 

Our method works equally well for random matrix theory. An interesting feature of this 
work is the new formula for the n-correlation of random matrix theory that arises by this 
method (see Theorem [51 ) It is a far less elegant formula than the usual determinantal 
expression, but it allows for direct comparison with the number theoretical result, illustrating 
the identical structure of the n-point correlations of Riemann zeros and random matrix 
eigenvalues. In fact, in the scaling limit (when the variables in the test function are multiplied 
by logT/27r and T, the height up the critical line, becomes large) then all of the arithmetic 
features of the formula for the n-correlation of the Riemann zeros disappear and it exactly 
matches our new formula for the the n-correlation of eigenvalues of unitary matrices in the 
equivalent limit. This identification allows us to prove that in the scaling limit the leading 
order terms for our n-correlation of the Riemann zeros have the expected determinantal 
form. See [TU] for the explicit derivation of the asymptotic limit in the case of the triple 
correlation of Riemann zeros. The higher correlations follow in exactly the same way. 

This point is significant in view of the difficulty in making this identification in other works 
on n-correlation and n-level density. In Rudnick and Sarnak this identification is proven in 
the case of test functions restricted to [—1, 1]; the proof is quite involved and makes serious 
use of the restriction on the support of the test function. Indeed this point forms a difficulty 
which shows up, for example, in the work of Gao [13] on n-level density for zeros of quadratic 
//-functions. Rubinstein [17] had evaluated this for test functions whose total support was 
contained in [—1, 1] and verified the determinantal form for functions restricted to this class, 
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analogous to what Rudnick and Sarnak did. Gao extended the range of support to [—2, 2] but 
for these test functions was unable to derive the determinant al form, due to combinatorial 
complexities. Here we handle the case in full generality without any mention of the test 
function. It is possible that our method will shed light on this difficulty that arises in these 
other works. 

This paper extends the calculation of the triple correlation of Riemann zeros [10J. An 
anticipated application of this current work is to the determination of the lower order terms 
in the nearest neighbor spacing for zeta-zeros. 

Throughout this paper we assume the truth of the Riemann Hypothesis. 

2. Background and notation 
2.1. The Riemann zeta-function. The Riemann zeta-function is defined by 

oo 1 
n=l 

for s = o + it with a > 1. It has a meromorphic continuation to the whole complex plane 
with its only singularity a simple pole at s — 1 with residue 1. It satisfies a functional 
equation which, in its symmetric form reads 

(2) ^ fr (i)^) = ^(nr)^ 1 - s) 

and in its asymmetric form is 

(3) COO = xWC(i - <0 

where 

(4) X (l - s) = X (s)- 1 = 2(2tt)-T( S ) cos H. 
The product formula discovered by Euler is 

p 

for a > 1 where the product is over the prime numbers p. 

The complex zeros of the Riemann zeta-function are denoted by p = f3 + ijj where it is 
known that < /3 < 1. The Riemann Hypothesis asserts that (3=1/2 for all zeros p. We 
assume this is true and denote the zeros as 1/2 + ijj, where < 71 < 7 2 < . . . . The number 
of 7 with < 7 < T is given by 

(6) N(T) = #{7 < T} = L log + O(logT) 

so that the average distance from one 7 to the next is ~ 2n / log 7. 

The family {£(1/2 + it)\t > 0} parametrized by real numbers t can be modeled by char- 
acteristic polynomials of unitary matrices. 



4 



J. B. CONREY AND N. C. SNAITH 



2.2. Unitary matrices. If X is an N x N matrix with complex entries X = (xjk), we let 
X* be its conjugate transpose, i.e. X* = where = a^-J- X is said to be unitary if 

XX* = I. We let U (N) denote the group of all N x N unitary matrices. This is a compact 
Lie group and has a Haar measure which allows us to do analysis. 

All of the eigenvalues of X 6 U(N) have absolute value 1; we write them as 

(7) e ie \e l62 ,...,e i(>N . 

The eigenvalues of X* are e~ l$1 , . . . , e~ ldN . Clearly, the determinant, det X = [j^=i e%6n °f a 
unitary matrix is a complex number with absolute value equal to 1. 

The average distance from one 9 to the next is 2ir/N. To obtain a sequence of numbers 
with average spacing 1 we let 

For any sequence of iV points on the unit circle there are matrices in U (N) with these points 
as eigenvalues. The collection of all matrices with the same set of eigenvalues constitutes a 
conjugacy class in U(N). Thus, the set of conjugacy classes can identified with the collection 
of sequences of iV points on the unit circle. 

We are interested in computing various statistics about these eigenvalues. Consequently, 
we identify all matrices in U(N) that have the same set of eigenvalues. Weyl's integration 
formula gives a simple way to perform averages over U (N) for functions / that are constant 
on conjugacy classes. Such functions are called 'class functions'. Weyl's formula asserts that 
for such an /, 

(9) / /(X) rfHaar = / f{B u 6 N )dX N , 

JU(N) J\0,2w] N 



where 



(io) dx N = 11 



e iVk _ e v 



,2 d9 1 ...dd 



N 



N\(2n) 



N ' 



Since N will be fixed in this paper, we will usually write dX in place of dX N . The Haar 
measure can be expressed in terms of the Vandermonde determinant 



(11) A(w u ...,w R )= det R (wj x ) = Y[ (w k - wj). 

The characteristic polynomial of a matrix X is denoted Ax(s) and is defined by 

N 

(12) A x (s) = det(J - sX*) = JJ(1 - se~ i9n ). 



n=l 
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The roots of Ax(s) are the eigenvalues of X and are on the unit circle. The characteristic 
polynomial satisfies the functional equation 

N N 

(13) A x (s) = (- a )"J] e -*»J](i- e *./ a ) 

n=l n=l 

(14) = (-l) N detX* s N A x *(l/s). 
Note that 

These characteristic polynomials have value distributions similar to that of the Riemann 
zeta-function and form the basis of Random Matrix models which predict behavior for the 
Riemann zeta-function based on what can be proven about the A. Some care has to be 
taken in making these comparisons because we are used to thinking about the zeta-function 
in a half-plane whereas the A are naturally studied in a circle. The translation is that the 
1/2-line corresponds to the unit circle; the half-plane to the right of the 1/2-line corresponds 
to the inside of the unit circle. Note that A x (0) = 1 is the analogue of lim^oo ((a + it) = 1. 
We let 

( 16 ) = 

In our formulas for averages of characteristic polynomials the function z(x) plays the role 
for random matrix theory that ((1 + x) plays in the theory of moments of the Riemann 
zeta-function. 

We want an accurate formula for 

(17) E* /(7*,---,7*J, 

0<7j 1 ,-,7j„<T 

for suitable functions /, to be described later, where the sum is for distinct indices j; the 
desired formula should be analogous to the RMT theorem which we state in the next section. 

3. Eigenvalue correlations 
Here is a statement for n-correlation of eigenvalues of random unitary matrices of size N: 
Theorem 1. Let f : [0,27r] n — > C be a continuous function of n-variables. Then 



I 



f( 9 ^ • • • > O jn )dX N = -i- / f(9 1 , ...,9 n ) det S N (9 k - Oj) d9 1 ... d9 n , 



u(N h< n ,..., jn <N ^ n > 
where indicates that the sum is for distinct indices and where 



sin ^ 
S N (9) = . % 
sin § 
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This theorem is a well-known consequence of Gaudin's Lemma, see 

In the following sections we present a new proof of this theorem, for periodic, holomorphic 
test functions /, based on a formula for averaging ratios of characteristic polynomials of 
unitary matrices. There are many proofs of this formula, see [H [3] . The point of this 
approach is that it has a natural analog in the theory of L-functions. 

3.1. Averages of ratios of characteristic polynomials. The statement of the ratios 
theorem is slightly complicated. We attempt to make it easier to comprehend by eliminating 
subscripts. So, let there be given finite sets A, B, C and D and consider 

(18) HA,B;C,D):= f ^i^^^ x, 

Ju(N) ll 7e c A M e 7 )lL e D A x*( e ) 

with 3?7 > 0, dtS > 0. Theorem [2], the ratios theorem, is written in an equivalent but slightly 
different form to previous work, where we express 7Z(A, B; C, D) as a sum over subsets 
S G A and T C B with \S\ = \T\. Each term in this sum essentially has the same structure, 
except that the elements of S effectively exchange places with those in T. In addition we 
let S = A — S and T = B — T. We will let a denote a generic member of S and (3 denote 
a generic member of T; we will use a and (3 for generic members of A and B or of S and 
T, according to the context. Also S~ = {—a : at G S}, and similarly for T". The Ratios 
Theorem is most easily stated in terms of 

(19) Z(A,B) := Y[z(a + (3), 



/3£B 



where z(x) = 1 _ e - x , and 



(20) Z(A,B;C,D) 



Yl«eAz( a + f3)U re czh + S) _ Z(A,B)Z(C,D) 
U^a z(a + 6) z(P + Z(A, D)Z(B, C) ' 

S£D 76C 



Theorem 2 (Ratios Theorem [TJ [8]). With > 0, > for 7 e C and 5 G D, \C\ < 

\A\ + N and \D\ < \B\ + N, we have 

TZ{A,B;C,D)= e~ N{j:& ^ &+j: ^P ) Z(S + T-,T + S~;C,D), 

SCA,TCB 
\S\ = \T\ 

where A = S + S , B = T + T and Z is defined at ( dZ| ). 

3.2. Averages of logarithmic derivatives of characteristic polynomials. For use in 

determining multiple correlation we differentiate the Ratios Theorem to obtain a theorem 
about averages of logarithmic derivatives: 

Theorem 3. IfStaij > and Ufa > for aj G A and /% G B, then J(A;B) = J*(A;B) 
where 

(21) J(A;B) := [ J] (-e^)^^) J] (-e^)^l(e^) dX, 

J W)ae A Ax Ax * 
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J*(A;B) 



(22) 



and 



E 

SCA,TCB 
\S\ = \T\ 



R 



Z*(S, S-)Z^{T, T~ 



(A-S) + (B-T) r=l 

=Ui+ ... +Ur 

\Ur\<2 



(23) H S , T (W) 



{a} C A 
{a, p} with 



S 
T 



imposing the additional 



E Ae 5T(«-«)-E^TT(« + / ? ) tfw 

Eteri(P-fi-Ea e si(P + &) tfW 
(i)'(a + (3) ifW 
otherwise. 

Also, Z(A,B) = TTaeA z(a + (3), with the dagger on Z\S,S~ 

/3GB 

restriction that a factor z(x) is omitted if its argument is zero. 

Remark 1. The definitions of J (A; B) and J* (A; B) make sense without the restriction that 
dlotj > 0, and ^t(3j > 0. However, the two are not equal without these restrictions. 

Remark 2. Note that J* (A; B) has a pole when an a G A is equal to —j3, for some (5 G B. 
It also appears to have a pole when two a's are equal, say a% = a<i, occurring when a.\ G S 
and «2 ^ S, as seen in the term —(a — a) of ( fffSj) . However, this is cancelled by a pole with 
residue of the opposite sign when S is replaced by S — {ai} + {c^} ■ The same phenomenon 
occurs when two (3 's are equal, as can be seen in the concrete examples given in $189\) and 



Proof. By (1181) . we have 



(24) 



J(A;B) 



C=A 
D—B 



Of course, in this situation \C\ = \A\ and \D\ = \B\; so that we may think of A = {a±, . . . , a^} 
and C = {71, . . . , 7^} and then the substitution "C = A" means the substitution 7, = for 
i = 1, 2, . . . , k, and similarly for D and B. 

Recall from Theorem [2] that 1Z is expressed as a sum of Z over subsets S and T. In 
performing the differentiations in (|24l) we will find that the derivatives with respect to the 
variables in S and T are fairly simple to perform (as we will show below, culminating in 
(I3ip ). but we will need Lemma [T] to differentiate with respect to the remaining variables. 
Hence we first rewrite Z so as to separate these variable types. 

Note first that Z(A, B) = Z(B, A) and that Z behaves nicely with respect to unions: 

(25) Z(A 1 + A 2 , B) = Z(A h B)Z(A 2 , B). 



Recall that A = S + S and B = T + T and put C = C$ + C-g and D = 
think of Cs, for example, as being the set that will be substituted by S when eventually C 



Dif where we 
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is substituted by A. Then using ( 1251) repeatedly, we have 
(26) Z(S + T~,T + S~;C,D) 

Z(S,T)Z(S,S-)Z(T-,T)Z(T-,S-)Z(C,D) 



Z(S, D)Z(T, C)Z(S~, Cs)Z(T~, D T )Z(S-, C S )Z(T-, D T ) 
This simplifies further if we make the substitution for C and D: 

Z(S,T)Z{S-,T- 



(27) Z(S + T-,T + S-;C,D) 



C=A 
D=B 



Z(S,T;T-,S-). 



Z(S,S-)Z(T, T~ 

Note that since z(x) has a pole at x = 0, the resulting expression is unless both S and T 
are empty. 

We now differentiate (|26|) with respect to the variables in S and T; these derivatives are 
easy to calculate because, anticipating the substitution of each 7 G Cs by an a we see that 
in differentiating with respect to a. the expression ,2(7 — 0:) in the denominator (one of the 
factors of Z(S~ , Cs)) must be differentiated; if not it makes the whole expression after the 
substitution is made because z(x) has a pole at x = 0. Using the notation 



(28) 

and noting that 
(29) 

we have, for example, 
(30) 

In this way we obtain 
(31) J(A;B) = 



Z\XX)= II z(x + y), 



x£X,y£Y 



d 



da ,2(7 — a) 



— e 



a— 7 



n 

&€S 



E 

S,T 
\S\ = \T\ 



d 



daZ(S-,C s ) 



c s =s 



(-1) 151 



n 

Per 



d d 



-N(J2 A+E P) Z(S,T)Z(S ,T ) 
Zt(5,5-)Zt(T, T-) 

Z(S, T)Z(5, S~)Z(T, T~)Z(C, D) 



da d(3 \Z(S, T)Z(Cg, S-)Z(D T , T~)Z(S, D)Z(T, C) 



C=A 
D — B 



Note that the sets C-§- and vary from term to term in the sum over S and T since the 
division of C into the union of sets Cs and C-g mimics the form of A = S + S, and similarly 
for D. Also observe that 



(32) 



Z(S, T)Z(S, S~)Z(T, T-)Z(C, D) 



Z(S, T)Z(Cs, S~)Z(D T , T-)Z(S, D)Z(T, C) 



1. 



C=A 
D — B 
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To perform the differentiations in (13T|) we use a form of logarithmic differentiation expressed 
in the following. 

Lemma 1. Let H be a differentiable function of w G W . Then 

(33) (Uj-) eH = eH E H( Wl )...H(W r ) 

\wew w J w=w 1 +-+w r 

where 

(34) H { W)=(t{{\h. 

The sum is over all set partitions ofW into disjoint sets Wj. 

In words this Lemma says that to perform a derivative with respect to each variable once, 
we form all of the set partitions of the complete set of variables and add up over these 
set partitions the product of the partial derivatives of the exponent H with respect to each 
variable in each subset of the partition. This lemma is obvious upon working a few examples. 

We apply this lemma with 

(35) H = H^°' D := J2 lo S *(« + P) + E lo S < a ~ d ) + E lo S *(P ~ P) 

aSS oSS /3£T 

/3GT SSS /jgT 

- J] log z(a + 5) - ^ log z(/3 + 7 ) 

QSS /3£T 
SeD 76C 

and so obtain, with H SjT (W) := H^ C ' D ( W) \ c =a = ((U w ew £) Hf^ C ' D ) \ c=a , 

' D=B \ ' / D=B 

(36) ^^Ee^'^ ^g'^, E f[»). 

•S| = |T| = 44^4 

Strictly speaking, ifs^iy) depends on A and -B, but from now on use of Hs,t(W) will refer 
always to the expressions in (|38|) - (H3|) and these can be used without specifically referring to 
A and B. 

By consideration of (135]) it is clear that we can restrict the subsets W r to be singletons or 
else pairs which have precisely one a and one (3. This follows from some easy calculations. 
Since 

(37) H^ c > D ({a}) = 7(0 + 0) + E ~ & ) ~ E + 6 ^ 

j3eT Sgd 



we have 

(38) H StT ({a}) = -(«- & ) ~ E ~(« + ^)» 



2 2^ 
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Similarly, 

(39) H s ,t({P}) = E^^-El^^)' P £T. 

/3eT «es 



In addition 



(40) H s ^ T ({a,(3})= I- J (« + /?), «,^5orT. 
Also, 

(41) # 5jT (0) = 1, 
and 

(42) #s,T({a, a'}) = #s,t({/?, /?'}) = 
and 

(43) H S , T (W) = 0, if |W| > 3. 



□ 



3.3. Residue identity. A key ingredient of the proof of n-correlation will be the following 
residue identity for J* (A; B): 

Lemma 2. Suppose that a* e A and f3* G B. Let A' = A- {a*} and B' = B - {/?*}. Then 
J* (A; B) has a simple pole at a* = —f3* with 

(44) Res J*{A-B) = NJ*{A'-B') + J*{A';B) + J*{A' + {-l3*};B'). 

a*=—/3* 

Proof. By Theorem [3] we have 

(45) J*(A,B)= £ D s , T (g,T) 

SCA 
TCB 

\S\ = \T\ 

where throughout this proof (and this paper) A = S + S, B = T + T, and 

(46) D StT (S,T) = Q(S,T) \[ H sAWr). 

S+T=J2W r r 

Here the sum is over any collection of non-empty sets Wi, W%, . . . which form a partition of 
S + T, 

(47) 0(S T)- e _JVC *« &+e ^t P) Z(S,T)Z(S-,T~) 

and Hs,t{W) is defined in ( I38l) - (l43l . We claim that D, Q and H have the following properties: 
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PI: If a* G S and (3* G T, then Q(S,T) is independent of a* and (3* and 

(48) Hsr{W) = {Et ,, + m 

v 1 ' v ' \ 0(1) otherwise 

P2: If a* G S and (3* G T, then Q{S,T) is regular when a* = -(3* and 

(49) ^W = {«^ + 0(1) T = { "* } 
v 7 ' v 7 \ 0(1) otherwise 

P3: If a* G S and /3* G T, then Q(S, T) is regular when a* = -f3* and 

(so) WstW = (^f + o(D if ^ = K> 

v 7 v 7 \ 0(1) otherwise 

P4: li a* E S and /3* G T and S' = S - {a*} and T' — T — {(3*}, then Q(S,T) = 
( R ^ F + 0(l))gi(-5,T) where 

Qi(S,T) = Q(S'X)(1 - (a* + f3*)(N + H s/ ,T'({a*})\ a * = -fi*+H s , :T ,({(3*})) 

(51) +0(|a* + /T| 2 )) 
and 

#S,tW = #S',T' W - («* + /?*) (#S',T'(W~ + K})a*=-/3* + ifs',T'(W~ + {0*})) 

(52) + 0(|a* + /T| 2 ). 

We show that the lemma follows from these four properties and then prove that these prop- 
erties hold in this situation. (We will later demonstrate a proof along very similar lines when 
we treat the n-correlation of the zeta-zeros.) 

From these four properties we obtain four Laurent or Taylor expansions of Ds,t(S,T) as 
a function of a* in a neighborhood of —f3*\ 

• If a* ES and (3* G T, then (with S' = S - {a*} and T' = T - {(3*}) 

DsAS,T) = ( ( ^ + 1 ^ )2 +Q(l))g(^T) U H sAW r ) 

v s'+T'=j:wr r 

(M) = ((^^ + 0(l)) flCT (S',r); 

consequently, Res a *=-p* D S ^ T (S,T) = 0. 

• If a* G Sand /3* G T, then' 

(54) ^Res^^^T) = Q(S' + {-/3*},T) £ J] 

= ^ ) 5'+{-/3*},t('S', T). 
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• If a* G S and p* G T, then 

(55) Res D SjT (S,T) = Q(S,T) V TT#S,T(W r ) 

s'+T=-£W r r 

= D s ,t(S',T). 

• Ha* ES and (3* G T, then 

(56) Qi{S,T) J2 i[H s , T (W r )=Q(S , ,T') ^ \[H s , jT ,{W r ) 

+ (a* + (3*) (N + H sl>T/ ({a*})\^ = ^ + H S ,, T ,{{(3*}) 

V #y,T'(^ + K})U» = -/3* + ffs',T'(Wr + {/?*}h 
^ #S',T'(W r ) / 

+0(\a* + (3*\ 2 )) . 



X 



Therefore in this final case, 
(57) Res D SjT (S,T) = ND s , jT ,(S,T) + D s , iT ,(S + {-/?*}, T) + Dg# r # (5, T+ {/?*}). 

a*=— p* 

Note that fl56l) can be written as 



(58) D s ,, T ,(S,T)(l + 0(\a? + 0*1)). 

By (1531) and (|56|) the double poles in PI and P4 cancel because 

(59) £ £> S ',T'(3,T)= J] D s ,T0,' i F) = r(A',B'); 

SCA,TCB SCA,TCB 
\S\=\T\ \S\ = \T\ 

{ a *}eS,{i3*}£T {a*}^S,{/3*}^T 

therefore, the pole at a* = — ft* is simple. 

Combining the four bullet-points above, we have (where as usual the primed notation 
means that a* or /3* has been removed from that set) 

Res r(A-B) 

a*=-j3" 



(60) = Yl (ND S >,T>(S,T)) + D s ,,t>(S + {-/?*}, T) + D s > iT /(S,T + {/?*}) 

TCS 

:|T| 

{/3*}£T 



SCA,TCB 
\S\ = \T\ 

{ a *}es,{0*}eT 



SCA,TCB SCA.TCB 
\S\ = \T\ \S\ = \T\ 
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Note that since a* doesn't appear in any of the summands, with the temporary convention 
that A' = R + R, we can relabel two of the sums as follows: 



SCA,TCB 
\S\ = \T\ 

{a*}es,{0*}er 



(61) 

Similarly 



Y D s , tT ,(S,T + {$*})+ Yl D s,t(S',T) 

tS SCA,TCB 

T\ \S\ = \T\ 

0*}eT { a *}£ S ,{/3*}eT 

= Y d rAW)+ Y d rAW) 

RCA',TCB 
\R\ = \T\ 

= J*(A';B). 



RCA',TCB 
\R\=\T\ 
{/3*}£T 



SCA,TCB 
\S\ = \T\ 

eT {a*}6S,{/3*}gT 

J*(A;B')\ a ^ = _^ — J*(A' + {—(3*};B'). 



SCA,TCB 
\S\ = \T\ 

{ a *}es,{p*}£T 



(62) 

Thus we have arrived at 

(63) Res J*{A;B) = NJ*{A'; B') + J*(A'; B) + J*{A' + {-/?*}; B'), 

a*=—f3* 

which is the statement of the lemma. 

Now we verify that properties PI through P4 are satisfied in the random matrix situation 
where we have 

Z(S,T)Z(S~ ,T~) 



(64) 
and 



Q(S,T) = e-^Eaesfi+EfcETfl 



Zi(S,S-)Zi(T,T~) 



{ E, eS 7(«-«)-E^7(« + ^) HW = {a}cA-S 
E^t E & es + HW = {P}CB-T 

otherwise 



if W = {a, (3} with °1*Zt 



(65) H s ,t(W) = < 

We will start with the simplest case that a* G S, (3* G T. The only polar term from 
a* = —(3* arises from a situation when one of the partition parts W r = {a*, (3*} and there is 
a pole from %(W r ) = (£)' (a* + (3*). Since (£)' (x) = 1/x 2 + 0(1) and a* and (3* don't 
appear in Q(S,T), this completes the proof of PI. 

Next, suppose that a* G S and (3* G T. The only pole in Ds,t(S, T) occurs in the product 
of the H for H s , T {W r ) when W r = {(3*}. We have 



(66) 



HsAin) = Y - - E -w + «) 

/5<GT <*es 
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for which, when a = a*, the term — ^{(3* + a*) has a simple pole at a* = —(3* with residue 
1. Q(S,T) is clearly regular at a* = —(3*. 

Next, when a* G S and f3* G T, the only pole in the product of the H occurs for 
H s , T {{a*})- We have 

(67) H s , T ({a*}) = £ £ + ^) 

* — ' z A — ' z 

for which, when ft — /3*, the term —^-(a* + (3*) has a simple pole at a* = — (3* with residue 
I. Q(S,T) does not depend on a*. 
Finally, we consider the case a* G S and /?* G T. We have 

(68) Q(S, T) = z(a* + P*)z(-a* - fl^Q^S, T) 
where 

Q 1 (S,T) = Q(S',T') 
m xe _jv (a . + ^) lifer' + - g) llagy f(6 + ~ 

Note that 

(70) Q 1 (5,T)| a , = _^=Q(5 / ,T'). 

We are well on our way to verifying property P4 because z(a* + f3*)z(— a* — (3*) = ^rp^pp + 
Y2 + 0(\a* + (3*\) and we have an expansion for Qx(S,T) in the neighborhood of a* = —(3*: 

Qi(S,T) = Q(S',T')(l-N(a*+[3*) + 0(\a* + f3*\ 2 )) 

x (l + (a* + /T)(E ( j(* + £*) - - &)) 

+ E (jW+^-J^-^))) + 0(|c** + /3*| 2 )) 
/3eT' 

(71) = g(5',T / )(l-(a* + /3*)(iV + i J5 , T ,(K})|^ = _ /3 , 

+f/^({ / 9*}))+0(|a* + /3*| 2 )V 
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Now we obtain an expansion for Hs,t(W), where we remember that W C S + T and so W 
does not contain a* or (3* . By (I65I) we have that 



(72) 



HsAW}) = ^-(ft-ft)-]T-(ft + /3) 



/3eT 



2 ^ 



(73) 

and 

(74) 

Thus, 

(75) 

and 

(76) 



/3er 



aes 



ft 



^5,T({«,/?}) 



^,t'({/3}) + -(/? - /?*) - + ft*); 
z z 



{a + p)=H 3 ,,T,({a,P}). 



H s ,t{W) 



H S ',t'(W) 



a*=-j3* 



d 

da* 



H S>T (W) 



-ft) (ft + /3*) if iy = {ft} 

-(7)' (/?-/?*) = {13} 

otherwise 



Note that we can write this as 
d 



(77) 



da* 



H S , T (W) 



a*=—(3* 



-H S ,, T ,(W + {a*})\ , - H S ,, T ,(W + {/?*}), 



where one or both of the terms will be zero. From this the expansion of Hs,t{W) in P4 
follows. 

This concludes the proof of Lemma [2j □ 



3.4. n-correlation via the ratios theorem. In this section we will prove the following 
expression for the n-correlation. 



Theorem 4. Let C_ denote the path from —5 + ni down to —5 — ni and let C+ denote the 
path from 5 — iri up to 5 + m and let f be a Ik -periodic, holomorphic function of n variables. 
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Using the notation J(A;B) from Theorem^ 

N 

/ Yl f(0n,---,0 jn )dX 
Jrx \n i» i 

(78) =-J- £ (-D'^'iV'-' 

^ ' K+L+M={l,...,n} 



x / / J(z K ; -z L )f(iz 1 , . . . ,iz n ) dz x . . . dz n 
Jc+ K JC- L + M 

where z^ = {z/~ : k G K}, —zi = {—Z£ : £ G L} and J c+K J c l+m means that we are 
integrating all of the variables in zk along the C + path and all of the variables in zl or zu 
along the C_ path. 

Proof. Since 

TV 

(79) g( z ) = K x (e z ) = X[(l-e z e-^) 

i=i 

has zeros at z = idj + 2mm, m G Z, by Cauchy's theorem we can express a sum 

(80) Y, W = T~- I - W/(^A) dz = h\ e z ^(e z )f(z/z) dz 
~ 2vr« Jc g 2m J c A x 

where C is a positively oriented contour which encloses a subinterval of the imaginary axis 
of length 2tt. We choose a specific path C to be the positively oriented rectangle that has 
vertices 5 — m, 5 + m, —5 + 7rz, —5 — 7ri where 5 is a small positive number. More generally, 
we have 

JV 

il,-Jn=l 



1 /" /" n A' 
™) Jc Jc fj[ A x 



v ' (2yri 

We average this equation over X G U(N) and, after a change of variables — > —Zj, we 
obtain 

r N 

(82) = — — — / J(«i,...,^;)/(izx,...,i^) d Zl ...dz n . 

(2m) n J C n 

Let C_ denote the path along the left side of C from —5 + 7ri down to —5 — m and let C + 
denote the path along the right side of C from 5 — m up to 5 + iri. Since the periodicity 
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of the function / implies that the horizontal segments of the contours cancel, each variable 
Zj is on one or the other of these two vertical paths. Thus, our expression is a sum of 2 n 
terms, each term being an n-fold integral with each integral on a vertical line segment either 
C_ or C + . For each variable Zj which is on C_ we use the functional equation ([15]) to replace 

A' A' 

e ~ Zj AA e ~ / by AT — e^^^-(e^). In this way we find that 
J (z u ...,z n ; )f{iz 1 , iz n ) dz 1 ... dz n 



(2m) 
(83) 



xf(iz±, ■ ■ ■ , i z n) dXdz\ . . . dz r< 
Another way to write this equation is 

/ J(z 1 ,...,z n ;)f(iz 1 ,...,iz n ) dz x ...dz r , 



(2m) 

(84) =7^1 (-1)" Ell/ ^^(e^)n/ U-e!»^(&) 
(2mYJ u{N) Kc ^ n} f e iJc + A x j£Jc-\ A*. 

xf(izi, . . . , iz n )dz\ . . . dz n dX. 

The expansion of the product over j ^ K can be easily expressed as a sum over further 
subsets of {1, . . . , n}. We have 

1 



(2m) 



J(zi, ...,z n ; )f(iz x , . . . , iz n ) dz x ... dz ri 



V ' J UyJN) K+L+M={l,...,n}jeK JL + A 0^ 



x II / N f( izi ' ■ ■ ■ ' iZn ) dzi. .. dz n dX. 

Using this last equation, fl82l) and the definition of J(A; B) from Theorem [3X we have the 
statement of the Theorem. 

□ 

3.5. n-correlation theorem. We will now prove our main theorem. 
Theorem 5. Let J* be as defined in Theorem^ Then 

I E f(9jiT--i9jn)dX N 

Ju ^ N h< n ,..., jn <N 

(86) = 77r^[ E NMj*(-i9 K ;i9 L )m,...,9 n ) dd x ... d9 n 



K+L+M={l,...,n} 
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where i9i = {i8e : £ G L}, —iBk = {—idk '■ k G K} and the star on the sum indicates sum- 
mation over distinct indices. Moreover, the integrand has no poles on the path of integration. 

Note the similar forms of Theorem H] and Theorem [51 In the former the sum is over all 
indices and the integrals are on paths slightly shifted away from the imaginary axis and in 
the latter the sum is over distinct indices and the integration is along the imaginary axis. 
Moving the integrals onto the imaginary axis results in some principal value terms, and 
surprisingly these cancel exactly with extra terms in the sum in [H 

We actually prove a more general theorem (Theorem [6] below). We start with a little 
notation: For a given n and < R < n, let the sum ■ , with the additional condition 

that j m 7^ ji if m > R and £ > R be denoted by ^ n,iJ . If we additionally fix three disjoint 
sets K, L and M whose union is {1,2,..., n}, then we introduce the following notation for 
the familiar integral 

/7ri pni p p 

■■ I I I J*{z K ; -z L )f(izi, iz n )dz x ■■■dz R dz R+l ...dz n 

-ni J-ni Jc^'-'V J c (L + M)n{l,...,R} 



17) = : I n > R 



+ 

f;K,L,M- 



Once again, the integrals on the imaginary axis are principal value integrals. 
We have already derived equation (1751) . In the new notation this is written as 



(2mT [ ^M:,..,^)^^ (-l) |i+M| ^ |M| ^,u- 

JU W K+L+M={l,...,n} 



Note that ([75]) features J whereas Ij'kl m ^ s defined in terms of J*. However, when R = n 
(that is, all the integrals are off the imaginary axis) Theorem [3] says that J and J* are equal. 
With the help of Lemma [2] we will prove the following: 

Theorem 6. Using the notation of (5?L) and the preceding paragraph, with < R < n, 



(89) {2m)» [ Y:' R f^---A n )dX= £ (-i)!^)^!....^}^^!^^^ 

JU ^ K+L+M={l,...,n} 



Proof. We will prove this by induction. Assume that Theorem [6] holds for n — 1 and any 
< R< n- 1. 

We start with the right side of ()89| and move the Zr integral onto the imaginary axis, 
resulting in a principal value integral and a residue at z R = z t , for t > R, in any term where 
R G K, t G L or R G L,t G K. A close inspection of the integral and the form of J*(zk', —zl) 
reveals that there is no pole unless R and t are in one of these two configurations (see the 
comment in the final paragraph of Section I37LI) . Also, if t < R then the contour on which zt 
is integrated has not yet been moved and so it remains on the far side of the imaginary axis 
from the zr contour and hence does not yield a pole. Each residue contribution comes in 
the form of the three terms in Lemma [21 multiplied by m. (It is iri rather than 2wi because 
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the zr contour is moving precisely onto the imaginary axis, where Zt lies, yielding half the 
contribution of a contour completely encircling the pole.) Thus 



(_iy(L+M)n{l,...,R}\ N \M\jn,R ^ ^ 

K+L+M={l,...,n} 

(90) = £ (_ 1 )l(L + M)n { l,.. > il-l}|^|M| 7 n,«-l^ 

K+L+M={l,...,n} 



n 

+2x m [ E (_ 1 )l(^+M)n{i,...,ii-i}| 7V |M| 

t=R-\-l K'+L'-\-M={l,...,n}-{R,t} 

+J*(z K >] -z L , +{t }) + NJ*(z K >; -^l')) 
x/(i^i, . . . , izR-uizt, iz R+1 , iz n )dzi ■ ■ ■ dz R _ x dz R+1 . . . dz n ^ 



The final sum above contains the two identical contributions from the case R G K,t G L and 
the case R G L,t G K. To confirm the sign of each term, if R G K, the residue is multiplied 
by +m because the contour of integration moves in from the right of the imaginary axis 
(skirting the pole in the positive direction) and the argument z R in J*(z K ; —Zl) occurs with 
a plus sign. Note that (_i)l(^'+M)n{i,...,R-i}| = (_ 1 )|(L+jion{i,...,fl}| if R e K and L = L > + { t }. 

On the other hand, if R G L then the z R contour comes from the left of the imaginary axis, 
but as C_ is directed downwards, the pole is still circled in the positive direction. However, z R 
appears in J*(zk', —zl) with a minus sign, so the residue acquires an extra minus sign, which 
is captured above because (_i)l(i'+Af)n{i,...,fl-i}| = (_i) x (-l)K i + M ) n { 1 --K}l if L = L' + {R}. 

In the integrals in the final sum above we now relabel the integration variables 
z±, Z2, ■ ■ ■ , zr-i, z R+ i, . . . , z n by zi, Z2, ■ ■ ■ , z n -i so that f{z\, . . . , z R -i, Zt, zr + i, . . . , zt, ■ ■ ■ , z n ) 
is replaced by f(zi, z R -i, z t -i, z R , . . . , z t -i, z n -i) =: g t {zi, z n -i). In addition, for 
some function h of sets K, L and M, 



(H K + {m}, L, M) + h(K, L + {m}, M) + h(K, L,M + {m})) 

K+L+M={l,...,m-1} 

(91) = Yl HK,L,M), 

K+L+M={l,...,m} 
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so we now rewrite the three J* terms in the final sum in ( |90l) as a sum over partitions of 
{1, . . . , n — 1}. Thus (jnOD equals 

E/_1 \|(L+M)n{l,...,ii-l}| aHM| rn.K-l 
V ' 1 f;K,L,M 

K+L+M={l,...,n} 

n 

(92) +2VTZ £ £ ( _ l) l(L + M)n { l > ...^l}|^|M| / n-^-l_ 

t=-R+l K+i+M={l,...,n-l} 

By the induction hypothesis, this equals 

El i N|(L+M)n{l,...,R-l}| a 7" I M | rn.fi-l 
V / 1 f;K,L,M 

K+L+M={l,...,n} 

n „ 

(93) +2tu E (27T0- 1 / E" :M 1 U.iO, 9 h MX. 

t=R+l JU W 

Note that the left side of ( 1891) can be written as 
(2m) 



JU(N) — 



(2ni) r 



(94) +(2mr / X) TT^StVh, ■ ■ ■ > 

where the second sum incorporates all the terms where 0.^ = 0j t , t > R, and then uses the 
same relabelling of the variables 9j 2 , . . . , 6^ _ 1? 6 , Jfl+1 , . . . , 6 , Jn and the definition of g t as 
described before ( |92|) . Therefore 

JU(N) — 

(95) = (-ifL+M^^R-l^^jn,^ 

K+L+M={l,...,n} 

and so, using the induction hypothesis, we have used (1891) for a given n and i? to deduce the 
same expression for n and R — 1. Since in (1551) we have derived the expression for R = n 
for any n, we have shown that if (1591 is true for n — 1, it is also true for n. To justify the 
induction hypothesis in n, we consider n — 1. Equation (175|) states 

(96) 2vu / X^/ftJdX = E (-l) l(i+M)n{1}l ^' M| 41^M = "^Jaw 



The final step above follows by remembering that J*($;za) = for any nonempty set A. 

-i,i _ _ 7-1,0 

/;0,0,{1} - J /;0,0,{1}' 



Since ' = X] 1 ' and il'l* (1) = — -^'00 m; h is immediate that 



07) 2ot 7 Y? fl w*) dX = E (-i) l(L+M)n0| ^ |M| ^^M- 

AM-L+A/={1} 
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This completes the proof of Theorem EJ □ 

It remains to verify that the integrand in Theorem [5] has no poles on the path of integration. 
We have already confirmed in Lemma [2] that each J* i#z,) has only a simple pole at 
k = -0 e for 9 k G 9 K and 6 e G 9 L . 

We check that 

(98) N M J*{-i9 K ;i9 L ) 

K+L+M={l,2,...,n} 

has no pole at 9\ = 9 2 for generic values of the remaining variables. A given J*(—iOK',i&L) 
only has a pole when 9\ G 9^ and 9 2 G 9k, or vice versa, so 

Res V N M J*(-i9 K ;i9 L ) 

01=02 ^— ' 

K+L+M={l,2,...,n} 

V N M Res (j*{-i9 K + {-i9 x }] {i9 2 } + i9 L ) 

K+L+M={3,...,n} 

(99) +J*(-i9 K + {-i9 2 }; {^x} + ^ L )) = 0; 



this is zero because Res s=x f(s,x) = —Res s=x f(x,s). Thus if (1981) had a singular set it 
would be of complex dimension less than n — 1 and this implies that there is no singular set 
(see for example [H], Corollary 7.3.2). 

Our new proof of n-correlation in the case of random matrix theory is now complete. 

4. Correlations of the Riemann zeros 

Now we turn to the Riemann zeta-f unction. The goal is to obtain a precise conjecture for 
the n-correlation of its zeros and we do this following the method of the previous section for 
the random matrix case. 

4.1. The ratios conjecture for the zeta-function. We derive our formula rigorously 
from the ratios conjecture for the zeta-function, which we now state. 

Conjecture 1 (Ratios Conjecture [6j). Let Z^(A,B) = [j^ ((1 + a + 0) and 

Further, let 

(101) A C (A,B;C,D) = T[Z P (A,B;C,D) [ A p o(A,B;C,D) d9 

P J o 

where z p (x) := (1 — p~ x )~ l , Z p (A, B) = Y[^a z p (1 + a + f3)~ l and 

(102) Z p( A,B;C,D)^ W)Z ' {C ' D) 



Z V (A,D)Z T (B,C) 
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and 

(103) AA A, B;C,D) := ^ ^"^f^ 

with Zp^(x) := (1 - e{9)p~ x )~ 1 . Then, provided that -\ < 3?a,3?/3 < \, < 3?7, < | 
an<i Sa, 9*/?, Q7, 95 -C T ; we conjecture that, with s = | + it, for any interval I C [-T, T] 



(104) 



a rue c(s + 7) lUd c(i - s + 5) 



ILeAC(s + a) rWC(l-* + l) 

h n 7ec c(«+7)n aeo c(i-«+5) 

w/iere 

(105) 7e fjt (A,S;C,D)= ^ T)%4 C (S + T~, T + S~; C, D). 



SCA,TCB 
\S\ = \T\ 



Here T means the set of all of the negatives of elements of T (i.e. T := {— t : t G T}), 
A = S + S, B = T + T and 

(106) X t (S, T) = J] X (s + &) J] x(l - s + $), 

&es p £T 

where %(1 ~ s ) — x( s ) 1 = 2(27r)~'T(s) cos ^ is t/ie factor from the functional equation 
C(s) = x(s)C(l-s). 

Remark 3. iVoie i/icrf since |5*| = \T\, for small shifts a and (3 we have 

(107) X t {S, T) = A + 0(1/(1 + j ? 

where i = log j- ; which can sometimes be used to simplify formulae. 

The method for constructing the ratios conjecture is detailed in [6] and is based on the 
same principles as the recipe for generating conjectures for moments (see [5]) of zeta and 
L-functions. (Moments cover just the case where C = {0} and D = {0}.) 

Corollary 1. With the same conditions on a,/3, 7 and 5 as in Conjecture [21 and with 
conditions on \i the same as those on a and (3, we have 

f EUa Cjs + a) Yhes C(l ~ s + (3) ^ x > 

Ji n 7£C as + 7) ca - * + 6) 11 / s + ^ 

(108) = / ft Cit (A, 5; C, D) ]J *-{s + fi) dt + 0(|/| 1/2+£ ) 
as a consequence of Conjecture^ 
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Proof. This follows immediately by integration by parts using the fact that 

(109) ^0) <log(2 + |s|) and ~^{a) < — Vi". 

X ds x !+ * 



□ 



4.2. Averages of logarithmic derivatives of the Riemann zeta function. To deter- 
mine the correlations of the Riemann zeros, we will need a result about averaging logarithmic 
derivatives of the zeta function: 

Theorem 7. Assuming the Ratios Conjecture, if$taa, 3fy3j > for ol% 6 A and f3j e B then 
J c ,i(A; B; U) = JIj(A; B; U) + O^ 1 ' 2 ^) where for an interval I 

(110) J U {A-B-U) 

■= fUj(i+ it + a ) n^i-^+^n (-^(i+^+m)) * 

and 

(111) 



J^; £?;[/) := J r Cjt (A;B;U) dt, 



Jl(A-B-U):= MS,T) 



SCA,TCB 
\S\ = \T\ 



z c (g,T)z c (g-,r-) 

Z C \S,S-)Z C \T,T-) 



A C (T-,S-;S,T) 



(112) 



x 



R , 

e n «s;t(w;) x n ( - -(V2 + « + a*)) • 



S+T r=l 
= Wl + -- + W R 



ifere we use the notation as in Conjecture^ and Z^(A,B) = Y[ C(l +ot + /3). In 

addition, T~ := {—t : t G T}, A = S + S , B = T + T and 

(113) H s , T (Wr) = H c , SiT {W r ) -Y, H p,i-MW r ) +J2H p , 2 . S)T (W r ). 



Further, we have 
(114) 



H ( ,s,t(W) = { 



Eaes C -(l + a-a)- J2per J i 1 + « + P) *fW = {a}c^ 
E deT C T (! + ~ P) ~ E & es T (1 + P + &) tfW = {(3} C T 



'per C 



i fW = {a,(3} with £g} 
otherwise; 
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H p ,i ;S ,t(W) = { 

and 
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Ea e *!(l + a-a)-Eter;*(l + « + /9) tfW={a}<zS_ 

EfeT t (1 + P - /?) - E & es J(l + P + «) ifW = {(3} C T 

f^Y (1 + or + /3) if W = {a, (3} with 

k otherwise] 



adS, 



(116) 



# P ,wW = J] (-l) J - 1 (J-l)!fJc S ,T(X,- 



117) c 5iT (X) := 



and i/ie notation 



(118) -V^T) :=^(T",5-;5,T), 

witt ^4 Pi 6i(v4, -B; C, D) as in Conjecture^ 

To prove this we want to differentiate TZ^j with respect to all of the a G A and f3 G B 
and then replace each 7 by an a and each 6 by a /3. In what follows A and C have the same 
cardinality, as do B and D. Thus, after differentiation, when we want to set the 7 equal 
to the a in some order, and the 6 equal to the (5 in some order, we can abbreviate this by 
C = A and D = B. With the definition of J<;,i(A] B; U) as in ffTTOl) . and using Corollary (TJ 
we have 



(119) J CJ (A;B;U) 



n 

/3es 



K c , t {A,B;C,D) 



n 



D—B fJ.£U 



- *-(l/2 + it + dt + 0{\I\ l / 2+e ) 
X ' 



The situation is much as in the random matrix theory case, except that now we have to 
understand how to include the arithmetical factor A(. For a start, we can differentiate with 
respect to the a G S and (3 G T as before. C$ and are defined as in the proof of Theorem 
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El We have 



(120) J C>I (A;B;U)= / JJ (- ^(1/2 + « + /*)) £ X t (S,T)^ 

" 1 aeU SCA,TCB C 

1 " |S| = |T| 



z c (s',T)z c (s , -,r-) 



(5,5-)z f '(r,r-) 



n 



/3ST 



gc(gz T )^' ^") Z c( T > T-)Z(C, D) _ 
L dadf3\ Z C (S, T)Z c (Cs, S-)Z C (D T , T~)Z ( (S, D)Z C (T, C) 



xA ( (S + T-,T + S-;C,D) 



C=A 
D — B 



dt + 0(\I\ 1/2+e ). 



In anticipation of applying Lemma [U we note that a brief calculation shows that 
Z d S, T)Z d S, S-)Z C (T, T-mC, D) + rJ + g _ ; ^ fl) 



D=B 



Z C {S, T)Z c {Cs, S-)Z C (D T , T-)Z C (S, D)Z C (T, C) 

(121) =^c(T- 5-;5,T). 

The remainder of the proof of Theorem [7] consists of applying Lemma [T] with (keeping just 
the factors from the big brackets in (I120p that depend on S and T) 

H = log Z C (S, T) + log Z ( (S,S~) + log Z ( (T,T-)- log Z C (S,D)- log Z C (T,C) 

(122) + ( T) + log 5") + log Z P (T, T~) - log Z p (S, D) - log Z P (T, C) 

v 

+ log I Ap, (S + T-,T + S-;C,D)d9. 
Jo 

Note the exponent —1 in the definition of Z P (A,B) in Conjecture [1] which accounts for the 
minus sign in front of H\ jPj s,T{W r ) in the definition in fl 113ft of Hs,r(W r ). 
Now it just remains to prove: 

Lemma 3. Let W C S + T. Then 

H p>2)S ,t{W) := H ^!og(^ A,e(S + T-,T + S-;C,D) dej 



c=s+s 

D-T+T 



(123) 

here 

(124) c s , T (X) :-- 



E (-i) J -V-i)!f[ c ^i); 



i=i 



d0 
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where we have adopted the notation 

(125) Ap t $(S,T) :=A pfi (T-,S-;S,T). 

Further, we have, for a* G S, (3* G T, S' = S - {a*}, T = T - {(3*} andW cS + T, 

d 

da* 

(126) +c s >,T'(W)c sl , T ,({a*})\ a * = ^ + c s , iTI (W)c s ^ T ,{{(3*}) 



-c s AW) 



= -c S 'AW + K})L*=-/?* - c s/ , T ,(W + {/?*}) 

-0* 



and 



<l H Pt2 . AT (W) 



da 



(127) -H p>2)S , jT ,{W + {a*})U*=-^ - H pAS , iT ,(W + {/?*}). 

Proof. The proof is simple, involving only differentiation. The first line, (I123p . follows from 
logarithmic differentiation of the integral of A Pi g, where each variable w G W appears in just 
one place. The equation (11261) also arises immediately by following the rules of differentiation. 
Note that a* appears in A Pt e{S, T) in the numerator of H103[) in a factor -2 P) 6»(| — a*) and in 
the denominator in a factor of type z p _q{^ + a*) and a* ^ W. The first term in (11261) comes 
from the z p _g(^ + a*) factor in the numerator of (I124p . the second term in (11261) from the 
z p,b{\ ~ a *) an d the final two terms in (11261) from the integral in the denominator of (I124p . 
To obtain (I127P note that the a* appears in each factor of cs,t(Xj), j = 1, . . . , J, in (I123p . 
Using the product rule on each term in (11231) we differentiate each c s ^r(Xj) in turn and sum 
the results. Note that cs^iXj) \ a *_ _** = cs',T>(Xj). Using (11261) and careful combinatorial 
accounting (11271) can be obtained. □ 

Remark 4. We have allowed here some loose use of notation in writing H P)2) s' \t'{W ' + {a*}) . 
The lemma starts out by defining H p ^ s> t(W) where W C S + T , but of course W + {a*} 
S + T when a* G S. However to understand the notation H p ^^s\T'iW + { a *}) simply replace 
S with S' and T with T' in the definition of H p ^,s,t o,nd replace S with A — S' and T with 
B — T' . A similar comment applies to cs',t'(W + {«*}) in U26]) . 

4.3. Residue identity revisited. Using Theorem U\ which is an application of Lemma [TJ 
along much the same lines as Theorem [31 we now prove the analogue for ( of Lemma [2J 

Lemma 4. Suppose that a* G A and (3* G B. Let A' = A — {a*} and B' = B — {(3*} and 
I = log t^. Then J£ t (A; B; U) (defined in Theorem^) has a simple pole at a* = —f3* with 

JEtes. Jl t (A; B;U)=-£{ 8 - (3*) J^A'; B'; U) 

(128) +r ( M'; B- U) + r Ctt (A' + B'; U). 
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Proof. We use Lemma |2j First we remember the convention that A = S + S and B = T + T 
and we write Theorem [7] as 

(129) Jl t {A,B;U) = H^{l/2 + it + v) ^ D C , S , T (S,T) 

*SJ X SCA 
\S\ = \T\ 

where, with the abbreviation A(S,T) := A$(T~, S,T), we have 

R 

(130) D CiStT (S, T) = Q C (S, T)A(S, T) £ J[ H s , T (W r ) 
and 

as!) Qds , T y., Xt{s , T) w^ri. 

Now we let Q((S, T)A(S, T) play the role of Q(S, T) in the proof of Lemma [2] and TCs,t(W) 
plays the role of H ST (W). Thus we need to prove the four conditions below, describing the 
behaviour of the various components of the formula as a* approaches —(3*, and then the rest 
of the proof is identical to that of Lemma [2j 

Ql: If a* G S and /3* G T, then Q C (S,T)A(S,T) is independent of a* and (3* and 

W ± w + 0{1) if W = R,/3*} 
0(1) otherwise 



(132) H s ,t(W) 



Q2: If a* G S and /3* G T, then Q C (S,T)A{S,T) is regular when a* = -(3* and 

v 7 v ' \ 0(1) otherwise 

Q3: If a* G 5 and /?* G T, then Q^(S,T)A(S,T) is regular when a* = -(3* and 

^V + 0(1) ifW = K} 



(134) H S , T (W) 



0(1) otherwise 



Q4: If a* G 5 and /?* G T and 5' = S - {a*} and V = T - {f3*}, then Q C (S,T) 
{j^Wr+0(l))Qt;i(S,T) and 



Q C;1 (S, T)^, T) = Q C (S', T')A(S', T) [l 

(135) -(a* + /?*)( - ^{s - (3*) + H s >,T>({a*})\ a * = -(3* + H S >A{P*}))) + 

X J 

and 

(136) H s ,t(W) = H s >,t>(W) - (a* + f3*){H S ',A w + K})U*=-/3* + H^W + {/?*})) 

+0(|a* + /T| 2 ). 
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We have just to prove these four conditions to complete the proof. We start with the 
first case where a* ^ S, [3* T. Then a* G S and (3* G T. The terms H P:1 and H p 2 have 
no poles because of the conditions on the real parts of a and (3. The only polar term from 
a* = —[3* arises from a situation when one of the partition parts is W r = {a*, [3*} and there 

is a pole from H C;S>T (W r ) = (1 + a* + (3*). Since (£)'(! + x) = 1/x 2 + 0(1) and 

Q^(S,T)A(S,T) is clearly independent of a* and (3*, the first condition is satisfied. 

Next, suppose that a* G S and (3* T. The only pole in D^. S ^ T (S, T) occurs in the 
product of the H for if^. SjT (W r ) when W r = {(3*}. We have 

(137) H c , s ,t{{(3*}) = £ £(l + - /T) - £ £(1 + /T + a) 

for which, when a = a*, the term — ^-(1 + (3* + a*) has a simple pole at a* = —f3* with 
residue 1. Q^(S,T)A(S,T) depends on a* and not (3*, so it is regular when a* = —(3*. 

Similarly, when a* ^ S and (3* G T, the only pole in the product of the H occurs for 
H C ,sA{®*})- We have 

(138) # C ;S,t(K}) = £ £(1 + & - a*) - £ £(1 + a* + /3) 

for which, when (3 = (3* , the term — ^(1 + a* + /?*) has a simple pole at a* = —(3* with 
residue 1. 

Finally, we consider the case a* e S and /3* G T. Let 5' = S - {a*} and T' — T — {(3*}. 
We have 

(139) Q c (5, T) = C(l + a* + /5*)C(1 - - Wc;i(S, 
where 

Q C;1 (S,T) = Q f (S',T')X t (K}, {/?*}) 

x lifer C(l + a* + 0KO- -<**-P) Yhes> CO- + & + glC(l -&-/?*) 
X IW C(l + a* - «)C(1 + a-a*) Ylp eT > C(l + - £)C(1 + 

Note that C(l + «* + /3*)C(1 - a* - /?*) = (a ,~^ )2 + 0(1). Also, remembering that 
X(s-/3*)x(l-s + /3*) = l, 

(141) Qc,i(S,T)\ =QdS',T>), 
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which gives us an expansion for Q^i(S,T) in the neighborhood of a* = —(3*\ 
Q Q1 (S, T) = Q c (S r , T')(l + ^(s - (3*){a* + (3*) + 0{\a* + (3*\ 2 ) 
x(l + (a* + {3*)(j2 (7(1 + « + /?*)- 7(1 -/?*-«)) 

+ E (7(1-^+^) -7(1 + ^* -/?))] + 0(|a* + /T| 2 )') 

/3eT' 

= Q C (S", T') (l - (a* +/?*)(- ^(s - /?*) + /W(K»U*=-/3* 
(142) +j ff C;S , T ,({^})) + o(|a* + /T| 2 )). 

Since ^.(S 1 , T)| a »__^» = A(S',T'), we have the expansion around a* = —(3*\ 

A(S, T) = A(S', T') (l + (a* + /?*) V ( V ( - ^(1 + a + /T) + ^(l-(3*- &)) 



+ E (--(i + 



0€T> 



+ 0(|a* + /3* 



= A{S', T') fl + (a* + (3*) (H PtliS ,,T>({<x*})\«-=-p + H^s'AiP*}) 
\ P 

(143) -H P!2 .s>, T >({a*})\ a * = ^* - H PtW ,({p*}))) + 0{\a* + (3*\)\ 

where the first line is a result of differentiating fllOip . and the second line from the definitions 
of H p>1 . SjT (W) (in TheoremED and H Pi2 . s>T (W) (in (jT251> ). 
Thus we have 

Q C;1 (S, T)A(S, T) = Q f (£', T')A(S', T 1 ) (l 

(144) -(a* + /3*)( - ^(s - (3*) + H s ,, T >({a*})\ a * = -p* + HvAiF}))) + 0(1) 

X J 
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Now we obtain an expansion for the product of H term. By the definition of in Theorem 
[7] we have that 



a — a) 



J2y(l + a + P) 



c 



(145) 



(146) 

and 

(147) 

Thus, 

(148) 

and 



(149) —H^riW) 



= #c;*',T'(M) + 4(1 + « - «*) - 4(1 + a + /?*); 

^C;5,t({/5}) = £ £(1 + (3 - $) - £ ^(1 + /? + «) 

/3er & es 

= F CiS , T ,({/3}) + £(1 + /3 - /3*) - £ (1 + + a*); 



H ( , S)T ({a,P})= (^j (l + a + P) = H c . t8 , >T ,({a,P}). 



H CiS ,, T ,(W) 



a*=-/3* 



(1 + a + P*) ifW = {a}cS 



a*=~l3* 



£) (1 + /?-/?*) HW = {p}cT 



C 

c 
c. 

otherwise 
-H C , S ,, T ,(W + {a*})\ a , = _p, - H CS >,T>(W + {/?*}). 



In exactly the same way, 

(150) H Ptl . StT (W) 

and 



H P)1;S >,t>(W) 



a*=-/3* 



d 



H p ,1:S,t(W)\ 



(151) = -H PiltSljTI (W + {a*})\ a *=^* -H pXS >,ti(W + {[3*}). 

Also, by Lemma [3] we have 



(152) 



H P ,2;S,t(W) 



H p ,2;S>,T>(W) 



a*=-/3* 
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and 

//,.2,sv/ 0n 



da* 



a*=-/3* 

(153) = -H P)2)3 ',T'(W + {a*})\^ = ^-H p , 2t s',T'{W + {ir}). 
Combining these results we have exactly equation (11361) . 

The rest of the proof proceeds exactly as before. □ 

4.4. n-correlation via the ratios conjecture. Now we proceed to n-correlation. Let / 
satisfy the conditions 

(154) f{xi, ■ ■ ■ , x n ) is holomorphic for \$sxj\ < 2, with j — 1, . . . , n, 

is translation invariant, ie. f(xi + t, . . . , x n + t) = f(x\, . . . , x n ) 
and satisfies /(0, X2, ■ ■ ■ , x n ) <C 1/(1 + \x2\ 2 + ■ ■ • + \x n \ 2 ) as \xj\ — > oo, 

with j = 2, ... ,n. 

Theorem 8. Let C_ denote the path from —5 + iT down to —5 — iT and let C + denote the 
path from 5 — iT up to 5 + iT and let f be as in ^15$ . Using the notation J^^(A\ B; C) from 
Theorem [7[ 

/(Tjh ■ ■ ■> Ijn) 

0<1 n ,-,J jn <T 

(155) = — V y (-i)i l i+i a/ i 

v ; (2ni ) n ^ K J 

v ' K+L+M={l,...,n} 

x / / ^ / J(A z k;-zl;-zm) dt f(izx, . . . ,iz n ) dz\... dz n 



c+ J 



K I p L+M 



where z K = : k G if}, — = {— ^ : £ G L} and j c K J c L+M means that we are 
integrating all of the variables in zk along the C + path and all of the variables in zl or zm 
along the C_ path; and I* is the interval which has lower endpoint max{0, — Qzi, . . . , —Qz n } 
and upper endpoint min{T, T — $szi, . . . , T — 

Proof. By Cauchy's theorem we can express the sum over zeros as 

/(7l>--->7n) 

0<7i,..., 7 „<T 

iff n C 

(156) = / • • • / /(-**!> ■ ■ ■ » - iz ») TT T^ 1 / 2 + Z A dzi '-- dZn > 

{zni) j c j c c, 

where C is a positively oriented contour which encloses a subinterval of the imaginary axis 
from zero to T. We choose a specific path C to be the positively oriented rectangle that has 
vertices 5,5 + iT, —5 + iT, —5 where 5 is a small positive number. 
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Due to the translation invar iance of /, (1156p equals 



, -IZn - t) 

c 



n 

X 

J 



"["J —(1/2 + Zj) dz\ . . . dz n dt 

1 ^ 



1 



f(-izx,...,-iz, 



X 



] J —(1/2 + it + Zj) d%\ . . . dz n dt 

j'=i ^ 



~1Zr, 



n „ 

(157) x JJ ^(1/2 + it + z,-) dt d2 t . . . dz n + 0(T e ) 

i=i ^ 

where the range of the innermost integral is the interval I* which has lower endpoint 
max{0, — ^Zi, . . . , — Qz n } and upper endpoint min{T, T — $Szi, . . . , T — Q^}. In the sec- 
ond line we made a change of variables Zj — > Zj + it. The contour C_ it is C shifted down by 
—it; that is, it runs from S — it, 5 + i(T — t), — 5 + i(T — t), — 5 — it. In progressing to the 
third line, we note that the horizontal portions of the contour of integration can be chosen 
so that the integral along them is 0(T e ) (following the identical argument to Davenport [TT] . 
page 108), so we concentrate on the vertical sides of the contours. When we now exchange 
the order of integration to move the t integral to the inside, the integration over z±, . . . , z n 
becomes the sum of 2 n integrals, each on one of the contours C + or C_ defined in Theorem 

E 

Remark 5. The main integral is of size ~ Tlog n T. The T is a result of j^_ TT ^ n f « T; the 
power of the log comes from the moment of the logarithmic derivative and will become clear 
from the examples at the end of the paper. 

For each variable Zj in (11571) which is on C_ we use the functional equation 

(158) C -(s) = ^(s) - - s) 

to replace ^(s + Zj), where s — 1/2 + it. In this way we find that (I157j) equals 

xf(izi, ■ ■ ■ , iz n ) dt dz\ . . . dz n . 
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Another way to write this equation is 

K ' Kc{l,...,n} j€K JC + j£K JC ~ J 1 S VA S 7 

xf(izi, . . . , zz n ) (it d^i . . . dz n . 

(Note that the dz's are no longer in order but this should not cause confusion.) The expansion 
of the product over j K can be easily expressed as a sum over subsets of K. This yields 

xjj/ ( - — (s + z m )) /(izi, • • • , dt dz 1 ... dz n . 
Remark 6. We note the asymptotic for 

(162) ^(1/2 + it) = - log K(i + (J^y 

In some applications \z m \ is small relative to t and it simplifies the formulae to replace 
^-(s + z m ) with — log^Jr. However, here where z k can be the same size as t we will not use 
this approximation. 

We have the statement of Theorem [SJ 

□ 

4.5. n-correlation for the Riemann zeros. We will now state our main theorem. 
Theorem 9. Assume the Ratios Conjecture^ Let J£ t be as defined in Theorem^ Then 

/(7i>---,7n) 

0<y 1 jt...jt 7n <T 



(163) 




xf( Zl ,...,z n ) dz x ... dz n + 0{T 1 ' 2 ^) 

where —izx = {—izk '■ k G K}, izi = {ize : £ G L}, and %zm = {iz m '■ m G M}. Moreover, 
the integrand has no poles on the path of integration. 

The proof is nearly identical to that of Theorem [5j The only difference is that some care is 
needed with regard to endpoints of intervals when we move each new path of integration onto 
the imaginary axis. The (slight) difficulty is with poles that may lie at the very endpoints; 
this point did not arise in the random matrix theory context because of the periodicity of 
the integrand. However by extending the paths slightly we can circumvent this difficulty; an 
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argument like that used to handle the horizontal segments in the proof of Theorem [8] will 
work in this case, too, and introduces an error term of size only 0(T e ). 

It remains to verify that the integrand in Theorem[9]has no poles on the path of integration. 
We have already confirmed in Lemma [2] that each J*(—i6x',iOL) has only a simple pole at 
9 k = -9 e for 9 k G 9 K and 9 e G 6 L . 

We check that 

(164) Yl Jtt(-^K;ie L ;ie M ) 

K+L+M={l,2,...,n} 

has no pole at 9\ = 92 when the values of the remaining 9j are unequal to 9\ or 9 2 . A given 
J^ t (—i9ic; i9l] i9m) only has a pole when 9\ G 9i and 6*2 G 9k, or vice versa, so 

Res V] j: t (-i9 K ;i9 L ;i9 M ) 

61=82 ' ' 

K+L+M={l,2,...,n} 

J2 fl Re f ( j U-Wk + {-i9i}; {i9 2 } + i9 L - i9 M ) 

' * 81=82 v 

K+L+M={3,...,n} 

(165) +Jtt(-^K + {-i0 2 }; {iOi} + i9 L - i9 M j) = 0; 

this is zero because Res s=x f(s, x) = — Res s=x f(x, s). 

Thus if (11641) had a singular set it would be of complex dimension less than n — 1 and 
by standard results in the theory of several complex variables, this implies that there is no 
singular set (see for example [Hj, Corollary 7.3.2). 

Our proof of n-correlation in the case of £-zeros is now complete. 

Corollary 2. By rearranging the integrals, now that we know the integrand has no singu- 
larities, and using the fact that f is translation invariant we have that the ratios conjecture 
implies that 

Q< 71 ^...^ 7n <T 

(166) = ^ f 77r~vi / Jl t {-iz K + it;iz L -it;iz M -it) 

1 J Q (2n) J[-t,t]* K+L+M={h n} 

Xf(zi, . . . , Z n ) dz X . . . dz n dt + 0(T l ' 2+t ) 

5. Examples 

In this section we explicitly write out all of the terms in our expressions for n-correlations 
for RMT eigenvalues and for £-zeros for 2 < n < 4. In the case of the ( correlations, we 
simplify the terms X t (S,T) and ^(s + a) using the approximations involving i = log ^ 
mentioned at ffTOTj) and (11621 . 

To proceed, we calculate J*(A; B) and J^ >t (A) B) := J£ t (A; B; U)/ ^-{\ + it + /i) for 

sets A and B with 4 or fewer elements. D$,t and D^ s,t are defined in the proofs of Lemma 
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|2] and Lemma HI respectively. In the following sections we first compile D$,t, Dc,s,t, then 
evaluate J* and J£ t and then assemble these into the correlation formulas, Rh )TI (xi, . . . ,x n ) 
and R^t,n(xi, . . . , x n ). Here we define R via 

/ V* f(9 1 ,...,9 n )dX N 

Ju ^h< n ,..,j n <N 

(167) = /0 \ T , / Rn,u(xi, . . . ,x n )f(xi, . . . ,x n )dxi ■ ■ -dx n 

{ 27T ) J[QM n 



and 



V fi.1l: • • • ) In) = To \n I ml R C,tA X l: 

(168) xf{x u x n )dx 1 ■■■dx n + 0(T 1//2+e ), 

(compare these with Theorems and ^ . 
Recall that 

(169) z(x) 



i •En) dt> 



Nr 

(170) S(x) = S N {x) - 



T 1 

sin | 



(171) z p {x) := {I - p-*Y\ 
and 

(172) Z P (A, B) = H z p (l + a + (3)~\ 

We also will introduce, as needed below, a number of other expressions A(x), B(x), etc.; 
for convenience, these are listed in Section 15.71 

5.1. Pair correlation, RMT. Suppose that the sets A and B have just one element: 
A = {a},B = {b}. We have 

(173) J*(a;b)=D^ + D atb , 
where 

(174) D M = Q (a + b), 
and 

(175) D af) = e- N(a+b) z(a + b)z{-a - b). 
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Thus, 

(176) J* (a; b) = Q (a + b) + e~ N ^z{a + b)z(-a - b), 
Then the 2-point correlation function is 

(177) R N ,2(u, v) = N 2 + J*(iu; -iv) + J*(-iu; iv) = det 



N S(u-t 
S(y -u) N 



5.2. Pair correlation, (. We have 

Jl t (a;b) = + a + 6) -£*(<* + &) 

(178) +e^ (a+fe) C(l + a + 6)C(1 -a- b)A{a + 6) 

where £ = log ^ and 



(179) 



(180) 



^ = 11 (Trip 



and conjecture that 

(181) Rc,t,2{u, v) = f + J^(iu; -iv) + J^(-iu; iv). 
Further, letting 

(182) P 1 (x) =e- ex A{x)({l + x)({l-x) 
and 



(183) P 2 (x) = (y)\l + x)-B(x), 
we have 

(184) J* t (a; 6) = Pi (a + 6) + P 2 (a + b). 

5.3. Triple correlation, RMT. In this case we have 

(185) J* (a; b u b 2 ) = D M + D aM + D aM , 

(186) = 0, 

(187) D aM = e- N ^z(a + h)z(-a - h) f-(b 2 - h) - -(b 2 + a) ) . 

V Z £ 
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and 

(188) D aM = e~ N ^z(a + b 2 )z(-a - b 2 ) (^-(h - b 2 ) - ^(61 + a)) , 
Thus, 

(189) J>; 61,62) = e~ N{a+b ' ) z(a + b 1 )z(-a-b 1 )(-(b 2 -b l )--(b 2 + a) 

\ z z 

+e -N(a + b 2 ) z{a + h)z{ _ a _ h) _ h) _ + a) 

z z 



Then 



Rn$(u,v,w) = N 3 + N(j*(iu; —iv) + J*{iv; —iu) + J*{iu; —iw) 

+J*(iw; —iu) + J*(iw; —iv) + J*(iv; —iw)) 
+ [J*(—iw; iu, iv) + J*(—iv; iw, iu) + J*(—iu; iw, iv) 

(190) +J*(iu; —iw, —iv) + J*{iv; —iw, —iu) + J*{iw; —iu, 

( N S(u-v) S(u-w) 
= det S(v -u) N S(v - w) 
\ S{w - u) S(w -v) N 

5.4. Triple correlation for (. The analogue for ( is 

J* t (a; 61, 6 2 ) = e-^ a+b ^A(a + 6i)C(l + a + &i)C(l - a - 61) 
x f ^(1 + 6 2 - 61) - ^(1 + a + 6 2 ) - Si (a + 61, a + 6 2 ) 

+e-^ a+b2 U(a + 6 2 )C(1 + a + 6 2 )C(1 - a - 6 2 ) 

(191) x ^(1 + 6 X - 6 2 ) - ^(1 + a + 61) - B 1 (a + 6 2 , a + 61) 



where 



+Q(a + bi, a + 6 2 ), 

log 3 p 



(192) QM = "d-^Xi-^) 



and 



(l-^)(l-^-jr + ^)logp 



IV 



38 J. B. CONREY AND N. C. SNAITH 

Then we conjecture that 

Ri;,t,3(u, v, w) = ^ + i(J^ t (iu; —iv) + J^ t (iv; —iu) + J^ t {iu; —iw) 

+J^ t (iw; —iu) + J^ t (iw; —iv) + J^ t (iv; —iw)) 

(193) + (j^ t (—iw;iu,iv) + J£ t (—iv;iw,iu) + J£ t (—iu;iw,: 

+J£ t (iu; —iw, —iv) + J£ t (iv; —iw, —iu) + J£ t (iw; —iu, —iv). 

It is convenient to introduce the function 

(194) P 3 (a, b, c) = B x {a + b, a + c) + ^(1 + a + c) - ^(1 + c - b). 

In terms of this, we have 
Jtt({a}; {h, b 2 }) = Q(a + h,a + b 2 ) - Pi(a + &i)P 3 (a, &i, h) -Pi(a + b 2 )P 3 (a, b 2 , h). 

5.5. Quadruple correlation, RMT. With A = {a},B = {b 1: b 2 ,b 3 }, we have that 
J* (a; 6i, b 2 , b 3 ) = D M + D aM + D aM + D aM 

= e- N ^z(a + b 1 )z(-a - b ± ) (^-{b 2 - h) - ^{b 2 + a)) (^-(b 3 - h) - + a)) 

(195) +e- N ^z(a + b 2 )z(-a - b 2 ) ( -(b, - b 2 ) - -{b x + a)) ( -(b 3 - b 2 ) - -(b 3 + a) 

\ z z J \ z z 

+e -N(a +b3 ) z{a + h)z{ _ a _ h) (^ {h _ h) _ *( bi + a) \ _ h) _ + a) 



With A = {a-i, a 2 }, B — {b 1: b 2 }, we have 

(196) J*(a>i, a 2 ] &i, b 2 ) = + D ai)bl + D aitb2 + D a2tbl + D a2)b2 + D aita2tbltb2 . 
where 

(197) D M = (^Pj + h) Q (a 2 + b 2 ) + {a, + b 2 ) (a 2 + h) 
and 

(198) D aiM = e - N ^ ai+bl) z(a 1 + b 1 )z(-a 1 -b 1 ) 

X (H {aih{bl} ({a 2 }, {b 2 }) + ^{a 1 },{6 1 }({a 2 })^{a 1 },{6 1 }({& 2 })) 

= e- N ^ +b ^z(a 1 + b 1 )z(-a 1 -b 1 )^^j (a 2 + b 2 ) 

+ ( -(a 2 - ai) - -(a 2 + 6i) J ( -(fo 2 - &i) - -(6 2 + <n) 
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the other D ai ^ j are similar. Also, 



f) _ p -A r (ai+a2+6i+62) v 

-"-^ai, 02,61,62 c 

z(ai + bi)z(—ai - b 1 )z(a 1 + b 2 )z(—a 1 - b 2 )z(a 2 + 6i)z(-a 2 - 61)2(02 + b 2 )z(-a 2 - b 2 ) 

z(ai - a 2 )z(a 2 - 01)2(61 - b 2 )z(b 2 - 61) 



Thus, 



(199) 



J*(ai, a 2 ; 61, 6 2 ) := - (ai + 61) - (a 2 + 6 2 ) + - (ai + 6 2 ) - (a 2 + 61) 



+e -iV(a 1+bl ) z ( ai + 6l ) z (_ ai _ &1 ) 



X 



- (a 2 + 6 2 ) + (-(a 2 
z I z 



+e 



-N(ai+b 2 ) 



z(a\ + b 2 )z(—ai - 6 2 ) 



I I (a 2 + 6i) + (^(a 2 
+e- JV(a2+6l) ^(a 2 + 6i)z(-a 2 - 61) 



+e 



x 

-jV(a 2 +6 2 ) 



- (ai + 6 2 ) + (-(ai 

2 / 2; 



^(a 2 + b 2 )z(-a 2 - b 2 ) 



x 



-J (ai + 61) + (-(ai 

2 / 2; 



ai) 



ai) 



a 2 ) 



a 2 ) 



-(a 2 + 6i))(-(6 2 -60 

z z 



-(02 + 62)) (-(61 -62) 

,2- ^ 



-(a 1 + 6 1 ))(-(6 2 -6 1 ) 

,2- ^ 



-(ai + 6 2 ))(-(6!-6 2 ) 

£ ,2- 



-(6 2 + ai)) 

z 



-(61 + al)) 

z 



— (b 2 + a 2 )) 

z 



-(6i + a 2 )) 

z 



z(-ai - 6i)z(-ai - b 2 )z(-a 2 - b 1 )z(—a 2 - b 2 ) 



x- 



z(ai - a 2 )2(a 2 - ai)z(&i - b 2 )z(b 2 - 61) 



40 J. B. CONREY AND N. C. SNAITH 

Then 



R Nj4 (u,v,w,y) = N 4 + N 2 (j*(iu;—iv) + J*(iv;—iu) + J*(iu;—iw) 

+J*(iw; —iu) + J*(iw; —iv) + J*(iv; —iw) 

(200) +J*(iy, ~iu) + J*(iu; -iy) + J*(iy; -iv) 

+J*(iv; -iy) + J*(iy; -iw) + J*(iw; -iy)) 

+N(j*(—iw; iu, iv) + J*(—iv; iw, iu) + J*(—iu; iw, iv) 

+J*(iu; —iw, —iv) + J*(iv; —iw, —iu) + J*(iw; —iu, —iv) 

+J*(—iw; iy, iv) + J*(—iv; iw, iy) + J*(—iy; iw, iv) 

+J*(iy; —iw, —iv) + J*(iv; —iw, —iy) + J*(iw; —iy, —iv) 

+J*(—iw; iu, iy) + J*(—iy; iw, iu) + J*(—iu; iw, iy) 

+J*(iu; —iw, —iy) + J*(iy; —iw, —iu) + J*(iw; —iu, —iy) 

+J*(—iy; iu, iv) + J*(—iv; iy, iu) + J*(—iu; iy, iv) 

+J*(iu; —iy, —iv) + J*(iv; —iy, —iu) + J*(iy; —iu, —iv)) 

+J*(—iy; iu, iv, iw) + J*(—iw; iy, iu, iv) + J*(—iv; iy, iu, iw) 

+J*(—iu; iy, iv, iw) + J*(iy; —iu, —iv, —iw) + J*(iu; —iy, —iv, —iw) 

+J*(iv; —iy, —iu, —iw) + J*(iw; —iy, —iu, —iv) 

+J*(iy, iu; —iv, —iw) + J*(iy, iv; —iu, —iw) + J*(iy, iw; —iu, —iv) 

+J*(iu, iv; —iy, —iw) + J*(iu, iw; —iy, —iv) + J*(iv, iw; —iy, —iu) 

( N S{u-v) S{u-w) S{u-y) \ 

S(v-u) N S(v-w) S(v-y) 

~ S(w-u) S(w-v) N S(w-y) 

\ S(y-u) S(y-v) S(y-w) N J 
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5.6. Quadruple correlation, (. We conjecture that 
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R(,ia( u , v i w , y) = ^ + ? 2 ( J it( iu ; ~ iv ) + J tt( iv ; ~ iu ) + J tt( iu ; ~ iw ) 

+J£ t (iw; —iu) + J£ t (iw; —iv) + J£ t (iv; —iw) 
(201) +JZ, t (iy, -*«) + Jlt(iu; -iy) + J£ t (ij/; -iv) 

+Jlt(iv; -iy) + J c * t (iy; -iw) + Jl t (iw; -iy)) 
+£(J£ t (—iw; iu, iv) + J£ t (—iv; iw, iu) + J£ t (—iu; iw, iv) 

+J£ t (iu; —iw, —iv) + J£ t (iv; —iw, —iu) + J£ t (iw; —iu, —iv) 
+Jl t (-iw; iy, iv) + J£ t (-iu; iw, iy) + J£ t (-iy; iw, iv) 

+ J tt(iy; - iw , ~ iv ) + J tt( iv 'i ~ iw i -w) + J tt( iw 'i -iy, ~ iv ) 

+Jl t (-iw; iu, iy) + J^ tt (-iy; iw, iu) + J c * t (-w; iw, iy) 
+Jlt(iu; -iw, -iy) + J£ t (iy; -iw, -iu) + J£ t (iw; -iu, -iy) 

+ J tt(-iy'i iu , iv ) + J it{- iv ' *y> iu ) + J it(- iu ; Wi iv ) 

+Jlt(iu; -iy, -iv) + J£ t (iu; -iy, -iu) + J£ t (ij/; -iu, -iv)) 
+J£ t (—iy; iu, iv, iw) + J£ t (—iw; iy, iu, iv) + J£ t (—iv; iy, iu, iw) 
+J£ t (—iu; iy, iv, iw) + J£ t (iy; —iu, —iv, —iw) + J£ t (iu; —iy, —iv, —iw) 
+J£ t (iv; —iy, —iu, —iw) + J£ t (iw; —iy, —iu, —iv) 
+J£ t (iy, iu; —iv, —iw) + J£ t (iy, iv; —iu, —iw) + J£ t (iy, iw; —iu, —iv) 
+J£ t (iu, iv; —iy, —iw) + J£ t (iu, iw; —iy, —iv) + J£ t (iv, iw; —iy, —iu) 

where the relevant J£ t are now described. 
We have 

(ono\ t* ( I n\ Ih h hl\ V" z p^ + a + h ^) z p( l + a + h *) z p( l + a + h) log 4 P 
l^U2) J ct ({a},{b 1 ,b 2 ,b 3 }) - - ^ 3+3a+bl+b2+b3 

v 

+W 1 (a, b x ; b 2 , b 3 ) + W x {a, b 2 ; b u b 3 ) + W^a, b 3 ; b u b 2 ) 

where 

(203) W^a, h; b 2 , b 3 ) = P^a + 6i)(P 3 (a, h, b 2 )P 3 {a, b u b 3 ) - B 2 {a, h; b 2 , b 3 )). 
with 

(204) B 2 (a,^M) = T (P-^'(P° +i '-l)(P° +i '-P)^P . 
v ' v ' ' ' ' ^ (-2p a+bl + p a+bl+1 + 1) (p bl - p b2+1 ) (p bl - p b * +1 ) 
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We also have 

Jlt({a u a 2 }, {h, b 2 }) = P 2 (a 1 + b 1 )P 2 (a 2 + b 2 ) + P 2 (a x + b 2 )P 2 (a 2 + 61) 

(205) -£4(01,02,61,62) 

+e _*( 01+02+ 6 1+62 ) A ^ ai ^ ^ y ^C({ fl l» fl 2>> (6l, 6 2 })^c({- a l ; ~ a 2}, {-61, -6 2 }) 

Z\{{ai, 02}, {-ai, -a 2 })^({6i, 6 2 }, {-61, -6 2 }) 
+W(ai, 61; a 2 , 6 2 ) + W(a u b 2 ; a 2 , 61) + W(a 2 , 6 X ; ai, 6 2 ) + W^(a 2 , 6 2 ; ai, 61) 

where 

(onrt\ A*tn „ h u\-TJ Z p(i a ^ Q 2)^ jK 6 2 })Z p ({-ai, -a 2 }, {-61, -6 2 » 

(206) A a 2 , 6 1; 6 2 ) - ^ ^ ^ ^ ^ _ M) 

Xp -a 1 -a 2 -6 1 -6 2 L + gpj 1 ~ a l ~ & l)^>(l ~ a 2 ~ 6l)z p (6 2 ~ &l) 



Zp(l)z p (-ai - bi)z p (-a 2 - h)z p (l + b 2 - bi) 
z p (l - ai - 6 2 )z p (l - a 2 - b 2 )z p (bi - 6 2 ) 
z p (l)z p (-ai - b 2 )z p (-a 2 - b 2 )z p (l + 61 - 6 2 ) 



and 

(207) W(ai, 61; a 2 , 6 2 ) = P^ + 6 1 )|p 2 (a 2 + 6 2 ) - B 3 (a u a 2 ; b u 6 2 ) 

+P 3 (ai, 61, 6 2 )P 3 (6i, ai, a 2 ) 

with 

(208) P 3 (ai,a 2 ;6 1 ,6 2 ) = ^log 2 p 

p 

+ 



Qr/ii _ p0 2 +l^ (_2p a i+ fo l + pOl+6i+l _)_ 1 ) 2 Qr>6i _ p 6 2+l^ 

C(ai,a 2 ;6i,6 2 ) 



(p<H _ ^a 2 +l^ ^_2p a i+ f >i -)- pai+61+l _|_ ^ (p 6 2+l — p b i) (pa2+&2+l _ 1) 
1 



pa 2 +f>2+i — 1 7 ' 



C( ai , a 2 ; 61, 62) = - P ai+bl + 2 P ai+bl+1 - p a ^+ 2 - p 2a 1+ 2fe 1+ l + p a 1+ a 2+ 2b 1+ l _ p a 1+ b 2 +2 
(209) _|_pa2+b2+2 _|_ ^201+61+62+1 _ 2p a i+ a 2+ 6 i+ 6 2+2 _|_ ^ai +02+61+62 +3. 

and 

^3 _pl+ai+6i _ pl+a 2 +6i _ pl+ai+b 2 _ pl+a 2 +b 2 _|_ p2+ai +02 +61+62 ) log 4 p 
P 4 (ai, a 2 , 61, 6 2 ) = (pl+ai+bi _ l)(pl+«2+6i) _ l)(pl+a!+6 2 _ l)(pl+a 2 +6 2 _ 1) ' 



CORRELATIONS OF EIGENVALUES AND RIEMANN ZEROS 



43 



5.7. Auxiliary functions. For ease of reference we list the various auxiliary functions we 
have introduced in this example section. 



(210) 



A ^ = 11 (Trip 



/ 011 x h h \ TT ^p({ fl i» fl 2}, {6i, 62»^p({-ai, -a 2 }, -&2> 

211 A tti.o^!,^ =11— — — - 

p z p(i a u a 2 }, {-ai, -a 2 })Z p ({bi, b 2 \, {-h, -b 2 }) 

Xp -ai -02-61-62 A + gpj 1 ~ Ql ~ ~ ° 2 ~ & l)*p( & 2 ~ &l) 



Zp(l)z p (-ai - b 1 )z p (-a 2 - h)z p (l + b 2 - b x ) 
z p (l - ai - 6 2 )%(1 - a 2 - b 2 )z p (bi - b 2 ) 



Zp(l)zp(-ai - b 2 )z p (-a 2 - b 2 )z p (l +h~b 2 ) J' 



(213) S 1 (x,?/) = 5] 



^ 1 ( 1 - ^ 1 1 1 - 2 + ^ 1p 2 -*+» 



P P 



ro14 , n( , , ,v (p - l)p2»i (p"+^ - 1) - p) log 2 p 

(214) B 2 (a,bi;b 2 ,b 3 ) = > = — — . 

V ' 1 ^ (_2p a+f>1 + 1) 2 (pbi _p6 2 +l) (p6i _p&3+l) 



(215) B 3 (a 1 ,a 2 ;b 1 ,b 2 ) = ^log 2 p 



(P- !) 2 (p ai+bl - l)V 1_ 



■6i 



(pai _ p<«2+l) (_2p a i+ b i + pai+6i+l _|_ I ) 2 (p6i _ p62+l) 



_^ C(a l7 a 2 ;b 1} b 2 



(pai _ pa2+l) (_2p a i+ fe i -)- pai+6i+l _)_ ]_) (p&2+l _p&i) (j/12+62+l _ X) 
1 



"^^02+62+1 _ 1 

(3 _ _ pl+a2+bi _ pl+ai+b-2 _ pl+a 2 +b 2 _|_ ^2+ai +02+61+62 ) log 4 p 

S 4 (ai, a 2 ; 61, 6 2 ) = 2^ (y+ai+6i _ l)(pl+a 2 +6i) _ l)(pl+ai+6 2 _ l)(pl+a 2 +6 2 _ 1) ' 



V 



C(ai, a 2 ; 61, 62) = -p ai+&1 + 2p ai+&1+1 - p a 2+6l+2 _ p 2ai+26i+l + p ai+a 2 +26i+l _ p a 1+ b 2 +2 
(216) _|_p a 2+&2+2 _|_ p2a 1 +b 1 +b 2 +l _ 2p a i+ a 2+bi+b 2 +2 pa 1 +a 2 +b 1 +b 2 +3. 



(217) Pi (2;) =e- ex A(x)((l + x)((l-x) 
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(218) P 2 (x)=(^)\l + x)-B{x), 

(219) P 3 (a, b, c) = Si(o + b, a + c) + ^(1 + a + c) - ^(1 + c - b). 



(220) Q(x,y) = 



log p 



(221) W(ai, 6i; a 2 , 62) = A(oi + &i){p 2 (a 2 + 62) - S 3 (ai, a 2 ; 61, 62) 

+P 3 (ai, 61, b 2 )P 3 (b 1 , ai, a 2 ) 



(222) Wi(a, 6 i; 6 2 , 63) = Pi (a + &i)(P 3 (a, 61, & 2 )P 3 (a, 61, 63) - B 2 (a, 61; fo 2 , 63)). 
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